Quark clusters in quark stars and possible 
astrophysical implications 
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Abstract 

A quark-cluster state, rather than the color-super-conductivity state, may appear in 
matter with low-temperature but high density, since the phase transition of chiral 
symmetry broken and that of color-confinement could not occur simultaneously. 
Such quark clusters might be stable in strange quark matter. Quark stars would 
have temperatures to be low enough to freeze by residual color interaction between 
the clusters, and become then solid ones. The charge numbers and the sizes of quark 
clusters, as well as the residual interactions in-between are investigated. The solid 
state properties of quark stars are constrained by observations. 
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1 Introduction 



Since E. Witten suggested that st range quark m atter may be the true ground 
state of color interaction system ()WittenL Il984r ) . it has been of great interest 
to find such matter in reality. There are hitherto many kinds of attempts, 
which can be classified into two ways. One is to detect directly strange quark 
matter on earth or in space (e.g., the ground and the AMS-II observations of 
strangelets in cosmic rays), the other is to study distant objects by means of 
astrophysics (e.g., to identify quark stars with strangeness). As for the nature 
of quark matter of strange stars, current researchers are focusing on quark 
condensation in momentum space (e.g., 2SC, CFL and LOFF states). Recently 
it has been proposed alternatively that quarks in strange quark matter with 
high density and low temperature may condense in position space (i.e., forming 
quark clusters which are not c olor-confined), a nd may be in a solid state when 
the temperature is much low (Xu, 20031 2004). This kind of solid quark stars 
could be applied to explain well the discrepancy between the observed glitches 



and free-precessions of pulsars (jXuL l2004t Zh ou et all 12004) . 
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In the following sections, based on the idea of quark clustering, observational 
constrains on the nature of quark matter with quark clusters are derived. 



2 Charge numbers and cluster sizes constrained by the featureless 
spectrum 



The 500ks Chandra observation of RXJ1856.5-3754 (|Drake et all 120021) shows 



a featureless Plank-like spectrum, which is not natural to be explained in the 
atmosphere models of normal neutron star. However, for the sim plicity nature 



the s pectrum might be fitted well in a solid quark star model (Zh ang et al 



2004). But if a cluster possesses more charges, it may behave like an atom 
which should have complex energy levels and would thus show correspond- 
ing atomic features in the spectrum. Because the observational low-band of 
thermal X-ray spectra is ~ 0.1 keV, this kind of "atomic features" can not 
be observed if the charge number per cluster, Z, is not greater than ~ 2. In 
order to explain the observational fact of featureless spectra, the number of 
quarks per cluster, n, can't be too large, because Z can be simply written 
as n ■ n e /(3n b ). In order to know the upper limit of n, we calculate electron 
density, n e , as well as baryon density rib near the surface of the star. As will 
be shown, these two parameters can be obtained b y solving the charge neu- 
trality equation and state equation of zero pressure ( Alcock et al. . 1986| ). and 



are determined by three physical constant from MIT bag model. 

The charge neutrality equation is 
2 11 

-n u - ~n d - -n s - n e = 0, (1) 

and the zero-temperature equation of state is 

^(n u + n d + n s ) -Y^(^i + Hin-i) = B, (2) 



where i (= u, d, s, e) denotes up, down, and strange quarks, and electron 
Qi, the thermal dynamic potentials; /ij, the chemical potentials ; and rij, the 



numb er densities. Here we take the formulae in the appendix of I Alcock et al 
(1986) to relate fi; and //j, 



2 



-^[3{„ - 8 {£ - m*yi* - ml H ^M^!1]Y - 2 {£ - m\f (5) 
_ m 4 l n 2 ^ + 6 In f^ml^l - m 2 ^ 2 - m\ l n [ ^±^±g^^ ]}]), 

n. - -i^- («) 



Beside eqs.(l)-(6), as a result of weak-interaction equilibriums between quarks 
and electrons, the chemical potential /i of them have the following relations, 

= Us = V, (7) 
ji u + fi e = (i. (8) 



We have also the relations between the number densities and the thermody- 
namic potentials, 

rii = - — . 9 
dm 



Using above equations, the number densities of quarks and electrons can be 
obtained, which are found to be dependent only on three parameters of MIT 
bag model: the vacuum energy density B, strange quark mass m s , and strong 
interaction coupling constant a s . Our results show that for some value of 
these parameters the electron density can drop to zero very fast. Beyond these 
parameter values, no solution for n e exists any more. Fig 1 shows the electron 
densities as parameters of m s and a s for different B. It is found that low m s 
and/or high a s favor low n e . This means that the electric field occurred on 
the quark-star surface could be very weak if m s is low and/or a s is high. 

With the charge density we now are able to constrain the number of quarks 
per cluster n from X-ray observation of RXJ1856. 5-3754. In order to reproduce 
the featureless thermal X-ray spectra, the maximum value of quark number 
in a clusters, n max , is calculated. Fig 2 shows the results of n max for different 
parameters in the MIT bag model (we choose Z = 2 in the computations). 
We find that the maximum quark number in clusters can vary in a very large 
range of values. For 50 < m s /MeV < 200 and < a s < 0.6, we have 10 4 < 
rw < 10 s . 
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Fig. 1. Contour plots of electron density in unit of cm 3 (solid lines) and of baryon 
number density in unit of fm -3 (dashed lines). The left panel is for i? 1 / 4 = 140 MeV, 
and the right for B 1 / 4 = 170 MeV. 




Fig. 2. The left panel shows the contour plot of n max for = 140 MeV, and the 
right for B 1 / 4 = 170 MeV. 

3 Cluster size and interaction strength limited by glitches 



It was pointed out bv IZhou et al.l ((20041 ) that if strange quark stars are in a 
solid state, pulsar glitches could be understood well by a sudden release of 
elastic energy accumulated when the shape of a star changes as its spin slows 
down. The elastic energy of a solid quark star is microscopically the sum of 
bond-energies of all links between clusters, and is determined by the strength 
of interactions and the sizes (i.e., n) of clusters. Unfortunately up to now the 
strength of interaction can't be calculated from first principles. We suppose 
thus its value to be of the same order of the chemical potential of quarks; that 
is within lMeV — lOOMeV. By input different strengthes of interactions by 
hand, we can obtain the magnitude of corresponding elastic energy. 



A strange star should be in a solid state before a glitch occurs due to increasing 
elastic energy. There are generally two kinds of movement that are against 
the solidification of strange stars; one is thermal motion, and the other is the 
quantum uncertainty (i.e., the zero-point energy). From the X-ray observation, 
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we know that the thermal energy may be at the magnitude of 0.1 — IKeV. The 
zero-point energy can be written as h 2 / (2ml 2 ) if quark clusters are moving non- 
relativistically, where / ~ [n/(3n(,)] 1//3 is the average distance between clusters. 
From these expressions we can find that, for small clusters which have less then 
thousands or hundreds of baryon numbers, the effect of quantum uncertainty 
could be larger than that of thermal motion. This zero-point energy should 
be much smaller than the depth (V p ) of potential well between clusters, V p 3> 
135n _5,/3 n^ 3 MeV, where is in unit of fm -3 and m ~ 300n MeV is suggested. 
We have V p ^> 5Qn~ 5 ^ 3 MeV if nb = 0.26fm~ 3 is chosen. This result means that, 
for n < 11, the bond-energy between clusters must be much higher than lMeV; 
otherwise the matter should become a quantum fluid. In an other word, if the 
interactions between clusters do exceed much from lMeV, then it is natural for 
quark matter to be solidified even when the clusters have only ~ 10 2 quarks. 



To ge nerate glitches, a solid quark star should have shear modulus ([Zhou et al 
2004) fi ~ 10 30-34 erg/cm 3 . The shear modulus fi can be derived from the 



Young's modulus E via fi = E/[2(l + u)], where v is called Poisson ratio, 
< v < 0.5. We see fi and E are of a same order of magnitude. For simplicity 
we can assume the potential of interaction is Lienard- Jones potential, V(r) = 
-A/^ + B/r^, and we have then / = (2B/A) 1 / 6 and V p = -V(l) = A 2 /(AB). 
For finding a linear coefficient which cause the Hook's law, we extract the 
Taylor expansion near the bottom of the potential well, 

K « = TW + W^ r - <T )1/6]2 + ° [r - (? T )1/6|S - (10) 



from which we find that the second order term can be written as 9/2 1 / 3 • 
A 7 ' 3 B~*' 3 (r — I) 2 . Using I and V p to eliminate A and B, the second order term 
becomes 36Vp(A///) 2 . At the same time, by considering a cylinder of length 
L and cross section S, the number of links between clusters at the direction 
of length in this volume is SL ■ 3rib/n; thus the value of elastic energy can be 
written as, 

E c ~^SLMV r (f?=\Y(^Y. (11) 



We can then easily get the Young's modulus Y = 72V p rib/n. The typical value 
of shear modulus is (1/3 ~ 1/2)Y in which ~ 0.26 fm~ 3 , consequently, 

fi ~ 10 34 rT^VMeV) (erg/cm 3 ). (12) 



As mentioned before, from observations, we have shear modulus fi ~ 10 30 34 
erg/cm 3 . The number of quarks per cluster, n, has resultantly an upper limit, 
which is 10 4 for V p = lMeV and 10 5 for V p = lOMeV. 
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4 Conclusions and Discussions 



We have obtained the upper limits of the quark number per cluster, n, from 
featureless spectra and pulsar glitches. The maximum value of n, n max , con- 
strained from featureless spectra is ~ 10 4-8 for parameters of 50 < m s /MeV < 
200, < a s < 0.6, and reasonable bag constants B. The upper limit of n for 
generating pulsar glitches is 10 5 if bond-energy between clusters is lOMeV. 
We conclude therefore that the upper limit of quark number per cluster could 
be ~ 10 4-5 in order for bare strange stars to reproduce both featureless X-ray 
spectra and pulsar glitches. 

What have been done in this paper are based on the solid quark star model. 
This model is only a phenomenological model, into which it is very necessary 
to research from first principles. Nevertheless, the conclusion of our paper may 
provide helpful information about the detail processes of quark clustering and 
matter solidifying. 

One problem in our paper is that we had used an equation of state derived from 
the traditional concept of free-quark matter, not clustered-quark matter. This 
problem may affect the result significantly if clustering changes the equation of 
state too much. Although this issue needs to be treated seriously, the deviation 
of equation of state may be slight, since color- confinement does not occur yet 
in the quark-clustering matter and the pr operty of the vacuum could thus be 



not much different from that described bv lAlcock et al 
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